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We investigate a f{R) modification of gravity tliat is exponentiai in tfie Ricci scaiar R to expiain 
cosmic acceleration. The steepness of this dependence provides extra freedom to satisfy solar system 
and other curvature regime constraints. With a parameter to alleviate the usual fine tuning of having 
the modification strengthen near the present, the total number of parameters is only one more than 
ACDM. The resulting class of solutions asymptotes to w = — 1 but has no cosmological constant. 
We calculate the dynamics in detail, examine the effect on the matter power spectrum, and provide 
a simple fitting form relating the two. 



I. INTRODUCTION 



Einstein's general relativity has proved to be a remark- 
ably robust theory of gravity. Large numbers of attempts 
have been made to modify it, e.g. for the purpose of ex- 
plaining cosmic acceleration without adding a separate 
dark energy component. Many of these have been found 
to have pathologies such as ghosts, unbounded energies, 
tachyons, or other ills [H . Gravity theories that are sound 
may still fail to provide acceleration, an early matter 
dominated period, or fall afoul of gravity constraints on 
solar system or other scales. 

Here we work within the framework of f{R) theo- 
ries that generalize the linear dependence on the Ricci 
scalar R in the Einstein-Hilbert action. Many such vi- 
able models exist but tend to have a restricted range 
obeying structure constraints, and in this regime are ef- 
fectively identical to ACDM as far as the expansion his- 
tory and distances. Furthermore these often take the 
characteristic curvature scale for the modification to the 
action to be that of the present matter density, i.e. when 
R = SnGptoday, so that the modifications kick in near 
today. 

We explore ways of ameliorating all three issues: eas- 
ing the restrictions of solar system tests and structure 
constraints on the curvature evolution, loosening the ex- 
pansion history, and relaxing the fine tuning. 

Considerable work has been done on f{R) theories; see 
0, Q for overall reviews and references. In this paper we 
follow most closely the setup of [1, H, @ . Our aim is to 
investigate the dynamics in more detail than usual, and 
relate it to the growth effects more explicitly, keeping 
close touch with the cosmological observations of both. 
In f|ll] we describe the model, its equations of motion, 
and the solutions for the expansion history and effective 
equation of state. The effects on systems with high cur- 
vature and density gradients, such as the solar system 
and galaxies, are discussed in i jllll along with the growth 
of structure. 



II. MODEL AND DYNAMICS 

The f{R) class of extensions to general relativity repre- 
sents a simple modification that preserves the maximum 
number of derivatives entering at two, and can have a 
well defined limit in which the theory reduces to general 
relativity, typically in the high scalar curvature regime. 
The gravitational part of the modified Einstein-Hilbert 
action is 
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where g is the determinant of the metric and R is the 
Ricci scalar curvature. 

In order to avoid including an implicit cosmological 
constant we require the low curvature limit of / to van- 
ish, i.e. /(O) = 0. At high curvature (e.g. in the early 
universe), we want the modification term to be negligi- 
ble compared to the usual term linear in R, so we require 
f{R » l)/i? ^ 0. We do not want to put in by hand 
a specific preference for dynamics at the present, i.e. a 
fine tuning, so we do not set the characteristic scale of 
/ in terms of Rq = SirGp^a = 1), where a is the cosmic 
scale factor and p the matter density. Rather, we give / a 
transition scale r and will fit for r given observations. Fi- 
nally, in order to satisfy the constraints from locally high 
curvature systems such as the solar system or galaxies 
that gravity must be very close to general relativity, we 
take / to have a steep dependence on scalar curvature, 
rapidly restoring Einstein gravity. 

This steepness will be a key ingredient in improving 
agreement with structure constraints while ameliorating 
fine tuning. Steep potentials have been considered pre- 
viously (see, e.g., [a, but here we explore in 
substantial detail the dynamics, from the deviation from 
matter domination, to the maximum equation of state 
variation, to the long time behavior in the equation of 
state function w(a) and phase space w-w' . We also give 
a quantitative analysis of effects on the matter power 
spectrum, as well as an accurate fitting form relating the 
expansion and growth effects. These results are of inter- 
est apart from the specific form adopted. We now seek 
the simplest form with fewest parameters yet freedom to 
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fit over a variety of characteristic curvature scales, not 
putting in special behavior just at the present density. 
All these conditions are satisfied by the ansatz 



variation, 



/(i?) = -cr(l-e-«/'-). 



(2) 



This simple form involves two quantities, and in fact 
we will see that the combination cr is equivalent to set- 
ting the dimensionless matter density today fim , so there 
is only one parameter besides 17^- Thus for this form 
the complexity is just one step more than the standard 
ACDM model. We focus on cosmological gravity from 
the matter domination through the accelerating eras, i.e. 
from high curvature to the asymptotic future; note pri- 
mordial nucleosynthesis will not be altered, as discussed 
in Sec. Ell 

The Friedmann equation of motion is modified to 0] 
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{H^Y + H^ 
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= rn^a ^ . 



(3) 



where a prime denotes d/d\na, a subscript R denotes a 
derivative with respect to R, H is the Hubble parameter, 
and p is the matter density, taken to be the only physical 
component of energy density. Note R = 12H^ + 3{H^Y , 
and for our exponential gravity model fjfn = {c/r)e~^'^ . 
Thus Jrr > 0, a critical stability condition (see [2| for a 
review) . 

Following with some modifications we define 
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(4) 
(5) 



to take out the leading order terms. The equation of 
motion ^ becomes two coupled first order equations 

x'h = xr/3, (6) 
1 XH + [d- 9(cmVr)a-*^]e-'"/'~ 
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l + a^{xH + d) m^fRR\6H^ 
where d = cr / {Qm?) and 
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11 + a^(4d + AxH + xr/3) 
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l + a^ixH + d) 



(7) 



(8) 



In this form we have carried out the explicit cancella- 
tion of several terms that nominally appeared dominant, 
making the numerical solution more robust. 

Treating the terms modifying the Friedmann equation 
([3]) as an effective dark energy density ilj;e(a), we can 
define the effective dark energy equation of state and its 
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In the future, we see that Qde ^ 1, i-e. the modifica- 
tions dominate and the effective dark energy density goes 
to the critical density. In the high curvature regime, / 
goes to a constant, —cr, and so one expects this to be 
related to f2A, or 1 — f2„i in the spatially fiat universe we 
consider. In general, the exact expression for cr is given 

by 
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so there is a defined relation between cr and flm- When 
the curvature is large, xh ~ e~^/^ (see below) and can be 
neglected and the product cr is explicit in terms of fl„i. 
Since cr corresponds to (all plots are for = 0.28), 
this leaves c as the only other parameter of the theory. 

The dynamics can also be understood fairly simply. 
In the high curvature limit when both xh and xr are 
small, w — > — 1, i.e. it starts from a frozen, cosmological 
constant state during the matter dominated expansion. 
Expanding the equations of motion at high redshift one 
finds that xh, xr ~ e'^/'' and are positive. So from 
Eq. (jlOp one sees that the field starts off phantom, i.e. 
with w < —1, and since x'^/xr ~ then from Eq. (|lip 
the field evolves quite rapidly, w' ~ a~^{l + w), much 
faster than a physical thawing scalar field which has w' — 
3(1 + U-) pill. 

Solving the equations of motion numerically, we dis- 
play the effective phase space dynamics of w'-w in Fig.[TJ 
We indeed see evolution from a frozen, cosmological con- 
stant state {w = —1, w' = 0) to the phantom regime, 
more swiftly than thawing ( "sublimation" ) , and then an 
oscillatory behavior where the field reaches a maximum 
value of |l-|-w|, then crosses back over w = —1, reaches a 
secondary maximum and quickly damps around w = —1. 
Note the second period of oscillations is highly damped, 
so small on the figure it appears as just a short dash 
around w = —1. 

In fact, there is an asymptotic behavior toward an ef- 
fective cosmological constant. The stationary point of 
the equation of motion is given by 



Xr = 0, 

XH - -d(c+l)(e«/'^ + c)-i 
with an implicit relation for R of 
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(13) 
(14) 
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So R settles asymptotically to a constant value, as does 
H, meaning that the future solution possesses w = —1. 
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FIG. 1: The phase space dynamics of the effective dark energy 
parameters w'-w is unhke any standard scalar field. It crosses 
between the phantom and normal parts of the phase space 
(often a signal of modified gravity), and "sublimates" rather 
than thaws from the initial frozen state. Successive oscilla- 
tions about A are barely visible here, showing rapid damping 
to the asymptotic cosmological constant. Thick portions of 
the curves show the evolution for a < 1. 



When the asymptotic value of e^^'" is much greater than 
one then R/m'^ 12d, or 30.86 for 17„ = 0.28. No vahd 
attractor solution exists for c < 1, because this would 
require R < 0. 

Figure [2] illustrates the numerical solution for R{a) for 
several values of c. We see that as promised R evolves 
very steeply due to the exponential form of f{R) and is in 
the high curvature limit e^/'' 3> 1 until recently. It then 
quickly reaches its asymptotic value given by Eq. . As 
c increases, r necessarily decreases to preserve cr, i.e. fim, 
and so e^^^ is always large and asymptotically R/m? = 
30.86. Thus, the evolution of R{a) is quite similar for 
all c ^ 3 (rjrr? ^ 5). The influence of the modification 
rises near the present for any allowed values of r, without 
additional fine tuning. 

To examine the expansion history in more detail, we 
consider the effective dark energy equation of state w(a) 
induced by the gravity modifications. Figure [3] shows 
this function for a variety of different values of c. As dis- 
cussed, the equation of state goes from w — —\ ■sX high 
redshift to w < —1 then crosses back to w > —1. The 
maximum departure from w — —\ decreases rapidly as c 
increases, going roughly as e"^"^ for large c. Thus as c gets 
very large, any deviations from ACDM are strongly sup- 
pressed. Figure [3] extends the plot to the future, showing 
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FIG. 2: The Ricci scalar curvature, normalized by the present 
matter density, evolves very rapidly from the high redshift, 
high curvature state to a constant asymptotic state (for c > 
1). As c gets larger than unity, R/m? approaches a constant 
given by 12 {Q,^ — 1) (i.e. 30.86 for Q,m = 0.28 as here) and 
the modification factor e^^^ always stays large. The influence 
of the modification occurs near the present for any allowed 
values of r, without additional fine tuning. 



the rapidly damped oscillation around w = —1. 

From the perspective of expansion history measure- 
ments, the distance to the cosmic microwave background 
last scattering surface agrees with the ACDM model with 
the same present matter density to 0.54% (0.2%) for all 
valid c (for c > 1.5). Distances to redshifts z < 2, e.g. as 
measured by the Type la supernovae magnitude-redshift 
relation, agree to 2.3% (1%) for all valid c (for c > 1.6). 
The parenthetical values for each correspond roughly to 
next generation observational limits. However, we will 
see in the next section that cosmic structure and its 
growth impose more severe limits on the parameter val- 
ues allowed. 



III. CURVATURE, STRUCTURE, AND 
GROWTH 

The f{R) class of models of most interest acts like gen- 
eral relativity at high curvatures, so as to preserve agree- 
ments with primordial nucleosynthesis and other early 
time physics. However they possess additional scalar 
degrees of freedom that can become apparent in lower 
curvature or high density gradient regimes [H]- Treat- 
ing as a scalar field one can define a Compton wave- 
length below which the effects on structure formation 



4 




FIG. 3: The effective dark energy equation of state evolution 
is shown for various values of c. As c gets large, the expansion 
history becomes indistinguishable from ACDM. 
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FIG. 4: The long term history of the effective dark energy 
equation of state is evolution from a cosmological constant 
state {w = —1), deviation to both the phantom (w < —1) 
and normal {w > —1) sides, and rapid damping to a future 
cosmological constant state. 



and bound structures become significantly different from 



general relativity. For example, the first parametrized 
post-Newtonian parameter, 7, takes the value 1/2 rather 
than unity as in general relativity. (In the high curvature 
regime inside the structure, general relativity is restored 
by the chameleon mechanism (isj . under conditions de- 
tailed by, e.g., 0, which the exponential model can sat- 
isfy. We address this further below.) 

The key indicators to these effects are the quantities 
fii and the effective Compton wavelength with respect 
to fn, usually written in terms of 0,0] 



B = 



1 + /fl {H^y ' 



(16) 



where B is the square of the wavelength in units of the 
Hubble scale. Note that 1 + > in the exponential 
model for all c giving a finite attractor solution, i.e. c > 1. 
This even holds during the radiation era since although 
we are used to thinking that i? = then, the value is in 
fact R <^ pi-ad - but one still has R ^ %-KGpm,oa^^ >• 
STrGpm.o ~ r so R/r ^ 1 at least back through pri- 
mordial nucleosynthesis and so 1 + /i?. 1. The con- 
dition B > 0, together with the vanishing of derivatives 
of / at high redshift, delivers a standard radiation era 
[J]. Given the scale £?, modifications to structure occur 
for wavenumbers k > aH/B^/^. This will begin to af- 
fect linear perturbation growth, k < 0.1/i/Mpc say, for 
B > 10~^. We will later calculate a more exact observa- 
tional constraint. 

The exponential form of Eq. ^ ameliorates the is- 
sues involved with deviations from general relativity in 
structure constraints, since fa and B ~ e~^^^. First, 
the evolution of the curvature and its effects are much 
more rapid than in many f{R) theories, so for a given 
present value of curvature, or fa, the deviation from gen- 
eral relativity can be much smaller in the past, when the 
structure formed, at a = 0.5 say. To rigorously quantify 
this argument requires numerical simulations of nonlin- 
ear structure formation in the specific model of modified 
gravity to robustly compute the "leaking" of the non- 
general relativity infiuence into the structure over time, 
as in [3| . If we can take advantage of this increased lati- 
tude, we can allow the characteristic curvature scale r to 
be more relaxed from the present matter density, giving 
greater freedom in the model. 

Figure [5] shows the steep dependence of fu and B with 
scale factor a and how this allows a range for the charac- 
teristic scale r, or equivalently the parameter c (related 
by Eq. [T^ . For example, in the n — 4 model of Q 
the value of B drops by a factor 3 from the present to 
a = 0.5 (redshift z = 1), while the exponential model 
achieves a drop of 300 (30,000) for c = 5 (c = 8). While 
the characteristic scale often seen in the literature is , 
e.g. /(i?) varies as (i?/m^)", here the dependence e"^/*" 
permits r to lie in the range r/m^ < 3.8 and still sat- 
isfy /i?(a = 0.5) < 10~^. Since is determined by the 
matter density today, this freedom somewhat eases the 
coincidence between the characteristic scale and today's 
curvature. 
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This steep dependence also somewhat eases the condi- 
tion for bound structures to recover general relativistic 
behavior: while the solar system, lying in the high cur- 
vature background of the galaxy, satisfies the conditions 
easily, our galaxy halo requires values of fu ^ 10^® (see 
0, Sec. HID for a detailed calculation). The sharp drop 
in fn when going to higher redshifts, as mentioned above, 
is stronger for the exponential model than most f{R) the- 
ories. Still, galaxy constraints will play a significant role 
and we should be cautious of values of c < 4 (imposing 
the condition at a = 0.5, or c < 7 if at a = 1) failing to 
live up to the chameleon mechanism. 

In f{R) models of modified gravity it is extremely diffi- 
cult to get an appreciable deviation from ACDM behavior 
on the expansion side of observations while still obeying 
constraints from the structure side. Using the very steep 
dependence of the modifications on curvature, and hence 
a rapid redshift dependence of the key structure param- 
eter fn (or B), exponential gravity can again do better 
than most models to ameliorate the situation. If we were 
to impose fn < 10~^, say, from structure constraints, 
this would imply a maximum deviation in the equation 
of state (occurring at the peak of the first phantom ex- 
cursion) |1 -I- wimax < 4 X 10~^ and c > 7 if evaluated at 
the present, but |1 -f- lulmax < 2 x 10"'^ and c > 4 if eval- 
uated at the time the structure forms, say a — 0.5. So 
even pushing all values as far as allowed, it is almost im- 
possible to get an appreciable signal in 1 + w. (Note that 
assuming the characteristic curvature modification scale 
is defined by the present matter density, i.e. rjrr? — 1, 
gives c — 15.4.) Interestingly, we find |1 + w|max is ex- 
cellently approximated by (l/2)i?(a ~ 1), relating the 
expansion and structure sides. 

We can check the magnitude of deviations from ACDM 
growth of structure, at least in the linear regime, by solv- 
ing for the evolution of linear density perturbations. We 
take into account both the change in the gravitational 
coupling strength (i.e. Newton's constant) and the scale- 
dependent effect of anisotropic stress in the /(i?) model 
(see, e.g., 0,[il]). Figure [B] shows the fractional deviation 
in the matter power spectrum as a function of wavenum- 
ber today. 

The deviations become noticeable on the //j Compton 
scale and have a magnitude determined by B. Due to 
the steepness of the exponential model, the deviations 
are not as severe for a given /ij(a = 1) as for shallower 
/(i?) models. In addition, the enhanced freedom in the 
characteristic scale allows models to be viable for values 
of c much lower than the c — 15.4 given by r/rn? = 1. For 
c > 5, the power spectrum deviations are less than 5% 
at fc = 0.1/i/Mpc, beyond which linear treatment must 
give way to numerical simulations. 

Interestingly, the steepness also quickly reduces the de- 
viations in the growth as one goes to higher redshift. The 
c = 3 model with the largest deviations shown (12% at 
the extreme), has less than 1% deviations in the power 
spectrum at a = 0.5. Recall that B can drop by two or 
three orders of magnitude in this model between a = 1 
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FIG. 5: The derivative fn (solid black curves) and the struc- 
ture parameter B (dashed red) , related to the Compton wave- 
length of the scalar part of the gravity modification, are plot- 
ted vs. c. The upper pair is evaluated at the present and 
the lower pair is at a = 0.5. Note that fa (and B) has a 
very steep dependence on scale factor a, so structures that 
might nominally have fn -ix 10~* today may have formed 
when /_R < 10"'' and so not exhibit observable deviations 
from general relativity. The maximum deviation of the effec- 
tive equation of state from a cosmological constant (dotted 
blue) closely follows B{a — l)/2 (dot-dash magenta). 



and a — 0.5. Note that because the expansion histories 
of c > 3 models are so close to ACDM, one can interpret 
Fig. [HI also as showing the deviation in growth from a 
general relativity model with the same expansion history 
as the f{R) models. We have verified this directly by 
turning off the gravity modifications in the source term 
of the perturbation equation: the maximum change is be- 
low 0.2%. Thus this figure shows the effects of modifying 
gravity for matched expansion. 

IV. CONCLUSIONS 

The idea of explaining cosmic acceleration without 
adding a separate energy density component is attrac- 
tive. Einstein's general relativity has proved to be highly 
resistant to modification, however, and in excellent ac- 
cord with observations. Theories generalizing the linear 
factor of the Ricci scalar curvature to a more general 
function, f{R) theories, are one of the main survivors for 
modifying gravity. 

In this paper we explored a simple model exponen- 
tial in the curvature. This has several beneficial conse- 
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FIG. 6: The matter power spectrum in f{R) theory deviates 
from the ACDM result in a scale dependent manner. In the 
exponential gravity case, a wide range of values of c have sup- 
pressed deviations and can be consistent with observations. 
The lowest curve has the same value c = 3 as the top curve, 
but plots the power spectrum deviation at a = 0.5, showing 
the effect of the steepness of the model. 



and qualitative results we present should be generally 
applicable to any sufficiently steep model however. 

The dynamics of the effective equation of state has 
several interesting properties, including "sublimation" 
rather than thawing from a cosmological constant state. 
Unfortunately, even the loosened bounds on from 
structure constraints do not allow the equation of state 
to have detectable deviations from a cosmological con- 
stant expansion history - one can view the failure for 
exponential models, steeper than any power law, as an 
indication of the difficulty in distinguishing viable f{R) 
theories' expansion from the cosmological constant. We 
also find an excellent fitting form for the maximum devi- 
ation from w = — 1 in terms of the Compton scale B, a 
direct relation between the expansion and structure char- 
acteristics. 

For structure growth, there exists an enlarged param- 
eter space that gives observationally allowed deviations 
from ACDM in the matter power spectrum. In partic- 
ular, because of the steep evolutionary dependence, fu 
today may be as large as ~ 10~^ and still agree with 
structure constraints. The sort of model discussed here 
is of interest in that we have increased freedom to have 
a viable f{R) modification of Einstein gravity that is yet 
distinct from standard ACDM. 
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